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SHEAE FLOWS IN MULTICELL SANDWICH SECTIONS 


By Stanley U. Benscoter 


SUMMARY 


Solutions are developed for shear flows in multicell sandwich 
sections for various cell- twist distributions . The problems of twisting 
and lending of a cantilever beam are first considered. More general 
cell-twist distributions are then considered, including the arbitrary 
distribution. A formula is also developed for the torsion constant of a 
sandwich section. 


INTRODUCTION 


The relation between load and deflection, for a homogeneous, elastic, 
isotropic plate, is expressed by a well-known fourth- order differential 
equation. A similar relation between load and deflection for a multicell 
sandwich plate should be given by a differential-difference equation 
since the internal shear flows and cell twists vaiy in discrete steps . 

In order to develop this differential- difference equation for the relation 
between load and deflection, it is necessary to know relations between 
bending moments, twisting moments, and deflections. It is convenient to 
replace the concept of twisting moment by the concept of cellular shear 
flows corresponding to the well-known shear flows in a multicell section 
according to the St. Tenant torsion theory. This paper is devoted to the 
determination of various relations between cellular shear flows and cell 
twists which may exist physically . Solutions are given for the distri- 
bution of shear flows required to preserve continuity of the warping 
displacements when the rate of twist varies from cell to cell in various 
manners . Symmetrical shear-flow distributions provide a resultant torque 
on the section but no resultant shear. Antisymmetrical distributions 
give neither a resultant shear nor a torque and, hence, are self- 
equilibrating. A formula is also derived for the torsion constant of a 
sandwich section with any number of cells . 


SYMBOIS 


A cell area 


c 


width of cell 
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d l mensionleBB factor defined "by equation (ih) 

Young's modulus 

shear modulus of elasticity 

height of cell 

moment of inertia of cross section about Z-axis 

torsion constant 

number of cells 

length of wall segment 

number of a cell (us ually an end cell) 

shear flow 

coordinate along cell wall 

wall thickness 

torgue 

displacement along cell wall 
axial (warping) displacement 
concentrated load 
number of an arbitrary cell 
vertical coordinate 
axial coordinate 

aspect ratio of wall segment (l/t) 
dimensionless factor defined by equation ( 17 ) 
dimensionless factor defined by equation (hi) 
proportionality constant for characteristic distributions 
Poisson's ratio 
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T shear stress 

(j> angle of twist 

6 rate of twist 

1 j f function defined by equations ( 7 6 ) and (JO) 


DIFFERENCE EQJATION RELATING SHEAR FLOWS AND CELL TWISTS 


A difference equation relating cellular shear flow to cellular twist 
(rate of twist) may be derived from, the condition that the warping dis- 
placement must be a continuous function around any closed path. This is 
the condition of continuity for thin-walled sections - Consider the sand- 
wich section shown in figure l(a) . The nth cell is Bhown in figure l(b) . 
The displaced position of the cell due to loads is shown by dashed lines . 
The angle of twist of the cell is indicated as 0 n . The width and height 
of the cell are indicated by c and h, respectively. The left web, AD, 
is web k and the right web, BC, is web k+1. 

The z-coordinate is taken parallel to the axis of a cell and Is 
considered to be positive in the direction extending outward from the 
paper. The s- coordinate indicates jthe distance along a wall seguent and 
is considered to be positive in the counterclockwise direction. Axial 
displacements (warping displacements) are indicated by w and transverse 
displacements along a wall segment, by u. These displacements are 
positive in the direction of positive coordinates . The shear stress in 
a wall sequent may be defined in terms of these displacements by 

The shear stress t, or the shear flow q, are positive when acting in 
the positive direction of the s-coordinate . The vertical web displace- 
ments y-jj. and yjj+i are considered positive upward. The angle of twist 
is given by 




gkKL - Fk 
c 


( 2 ) 


The rate of twist is obtained by differentiating equation (2) with 
respect to z. Thus, 




( 3 ) 


In the top and bottom walls the quantity dufdz may be considered to 
consist of two parts. The first part is due to the rotation of the cell 
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and is given by the formula 0 n h /2 • The second part is due to the 
development of tensile or compressive strains and may he indicated as 
u-^(s) . This function is antisymmetrical about the center of the wall 

segment. In the vertical webs the displacement u is equal in magnitude 
to the displacement y. In the right web, u equals y-^-^; and in the 

left web, u equals -y-^. 


Equation (l) may be written 


for each web separately as follows: 


Web AB' 


Web BC 


Web CD 


Web DA 


T = G If + 6n + U1 ^ 


t = G |a + G 

ds oz 


T 



Gh a 4. 
~2 e n + 


ui(s) 


T 


g|2 - 

os oz 


(4a) 

(4b) 


(4c) 

(4a) 


The condition of continuity of warping displacements is expressed by the 
following equation: 


* 



(5) 


The subscript on the integral sign indicates that the line integration 
is carried around the nth cell. The equation must hold true for each cell. 


Integration of equation (l) around the nth cell gives, with substitu- 
tion of equation ( 5 ) / 




( 6 ) 


From the formulas in equations (4) it is seen that the integral on the 
right may be evaluated to obtain 


(7) 
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The function u^(s) makes no contribution to the integral . From. 

equations (2) and ( 3 ) the second term on the right may be evaluated to 
obtain 

^ n T ds = 2AG0 n (8) 

Introduction of the shear flow q = tT gives 

fn. ^ f ’ 

The total shear-flow distribution may be considered to consist of 
two parts . The first part consists of shear stresses associated with 
bending action. These shear stresses are related by equilibrium conditions 
to the rate of change of normal stress . The second part of the Bhear-flow 
distribution consists of cellular shear flows resembling the shear flows 
of multicell sections in elementary torsion- The shear flow in the top 
or bottom wall of a cell is equal to the cellular shear flow associated 
with the cell- The shear flow in a vertical web is equal to the 
difference between the cellular shear flows associated with the two 
adjacent cells - 

In the- present paper only the second part of the shear-flow distri- 
bution is considered. The shear flows accompanying bending are not 
considered further in order that attention may be devoted to a development 
of formulas for cellular shear flows with variable cell twists. 

Equation (9) becomes 


-°2Vl + ^ a)q n _a 2Vl = 


( 10 ) 


where 



( 11 ) 


The coefficient Sx is the sum of the "aspect ratios" of the wall segnents 
of a single cell. Structural dimensions of several sandwich sections are 
illustrated in figure 2. Although the equations ' have been developed for 
a rectangular cell, they are also applicable to other shapes. The 
structures of figures 2(b) and 2(c) are mathematical! y Identical insofar 
as the relationship of equation (10) is concerned- 


In the development of solutions of difference equations, it is 
generally convenient to use the letter x for the independent variable. 
Thus the arbitrary, cell is indicated as cell 2 . The shear flow of this 
cell is q 2 and the cell twist is 0 X - The last cell may be indicated 

as cell n . The number ing of cells is shown in figure 3 • A symmetrical 
numbering method may be used as shown in figures 3( a ) and 3(b) - 
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For an odd number of cells the total number of cells is 2n + 1, -while 
for an even dumber the total number is 2n. With either an even or odd 
number of cells the numbering from left to right as shown in figure 3(c) 
may be used. In this case, h is equal to the number of cells. With x 
as the independent variable, equation (10) becomes 

-^Vl + " °2Vl * (12) 


U0TFORM DISTRIBUTION OF 0 (PURE TORSION) 


Consider a cantilever beam as Bhown in figure k loaded with a 
torque T at the tip. According to the St. Venant theory all cells of 
an intermediate section experience the same rate of twiBt. The cell 
twist 9 X is replaced by a constant twist 0, and equation (l2) becomes 

” a 2 q x-l + a 2 q xfi 0 2Aae (!3) 

It Is frequently found convenient in expressing the solution of a system 
of linear algebraic equations to introduce a symbol for the ratio of 
nondiagonal to diagonal coefficients. In the present case the following 
factor is useful: 

' ( 14 ) 

Methods for solving difference equations of the preceding type are given 
in references 1 and 2. 

The solution of equation (13) may be expressed in terms of exponential 
or hyperbolic functions . The latter prove to be most convenient in the 
present case. The cells are considered to be numbered in a symmetrical 
manner as shown in figure 3(a) • The complete solution of equation (13) 
may be expressed in the following foxm: 

qx = B sinh f3x + C cosh 0x + D (15) 


where 


D = 


2AG0 


(1 - 2d)£a 

and the section property {3 is defined by the equation 


( 


(16) 


( 17 ) 


This solution may be verified by direct substitution into the equation. 
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The first two terms form, the general solution of the homogeneous equation. 
The last term is the particular integral. 

The constants of integration, B and C, must he determined. Because 
of symmetry the coefficient of the antisymmetrical term must he set equal 
to zero. Thus, 


B = 0 


and 


q^ = C cosh Px + D (iB) 

In order to determine the second constant of integration, it is necessary 
to define .a boundary condition. This boundary condition may he determined 
by considering the end cell (nth cell) of a section as shown in figure 5 . 
The shear flow in the top and bottom walls is q^ A comer element of 

the wall has been illustrated by cross-hatching at the . top and bottom of 
the end wall. In order that these corner elements shall be in equilibrium 
Jn the axial direction, the shear flow in the end wall must also have the 
value q n - An imaginary cell may now be considered as shown by the 

dashed lines in figure 5* This cell would be numbered n + 1. In order 
that the comer elements remain in equilibrium, the shear flow in the 
imaginary cell must be equal to zero. This may also be conveniently 
stated as the condition for preserving continuity of shear flow around 
the exterior comer of the cross section. The assumption of zero Bhear 
flow in the imaginary cell corresponds to the zero value of Prandtl's 
stress function at the edge of a solid member in torsion. The boundary 
condition is given by the equation 

<3 n +l = 0 (19) 

This boundary condition was made known to the author by M. A- Biot of 
Brown University. The constant of integration is determined by writing 
equation ( 18 ) for the imaginary cell. Thus, 

C cosh P (n + l) + D = 0 


or 


C = 


-D 


cosh p (n + l) 


( 20 ) 


Substitution of equation (20) in equation (l8) gives the solution 

cosh Px 


= D 1 - 


coBh P (n + l) 


(x = 0 , 1 , 2 , ...) (2D 
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When there is an even number of cells, they must he numbered as shown 
in figure 3 ("b) • This peculiar numbering is necessary in order to 
preserve the finite difference of unity between adjacent cells thro ugh 
the plane of symmetry • The boundary condition is that q j. be equal 


to zero • The resulting solution is 


= d 


1 - 


cosh Px 


cosh p^n t 


c X = \> \> \> ( 22 ) 


A graph of the shear flows is shown in figure It (c) for a seven- 
cell member. The graph may be regarded as a section through the Prandtl 
stress surface. The ord i nates to the steps, or plateaus of the stress 
surface, are the cellular shear flows. The cells have been considered 
to be square with the same thickness for all wall segments. Thus, 


and 


and consequently 


ti = tg 


d = 0.25 
cosh p = 2 

P = 1.317 


The magnitudes of these shear flows are as follows: 


^0 

q l 

q 2 

q 3 

0 .990D 

0-979D 

0.928D 

0.732D 


LUJEAR DISTRIBUTION OF 9 (BEAM BENDING) 


Consider the cantilever beam of figure 6(a) as being acted upon by 
a tip load W. It has been shown by Goodier (reference 3) that the 
St. Tenant bending theory leads to the conclusion that cell, twists at an 
intermediate section have an anti symmetrical linear distribution. The 
corresponding shear flows form a completely self-equilibrating system. 

The cell twists that arise are associated with the pronounced anticlastic 
shape of the deflection surface which occurs in wide beams of shallow depth . 
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The cell twist may be expressed as a linear function of 2, that is. 


e z = -kx 


(23) 


Keference 3 shows that the constant k, for team tending, can he expressed 
as 


k = 


pWI, 


EX 


S 


(2k) 


The values of Poisson's ratio p and Young's modulus E are assumed to 
he the same for all wall elements . The moment of inertia of the total 
team cross section is indicated as I Z' 


The difference equation (equation (l2)) becomes 

-a 2 q x-l + " °2Vl = _2AGlac (25) 

The general solution Is given ty 


q 2 = B slnh Px + C cosh 0x + Ex 


where 


F _ -2AGk 
(l - 2d)Zo. 


( 26 ) 

(27) 


and P is defined ty equation (17) • Since the solution must be anti- 
symmetrical. 


C = 0 


and 


q x = B sinh, Px + Ex 


(28) 


The boundary condition must again be a vanishing of the shear flow in an 
imaginary exterior end cell. Thus, 


^n+1 = 0 


(29) 


Therefore, equation (26) becomes 


q z = E 


x 


(n + l)sinh Px 
sinh p(n + 1) 


(x = 0, 1, 2, ...) ( 30 ) 


A solution of this type was given in reference 3 but differed somewhat 
because of the different boundary condition in the example considered. 
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A graph of these shear flows is shown in figure 6(c) . The magnitudes 
of the shear flows for a seven-cell section, with square cells and constant 
wall thickness, are as follows: 


q o 

q l 

q 2 

q 3 

0 

0 . 928 F 

1.7I1F 

1-918F 


SINE DISTRIBUTION OF 0 


Several solutions will now "be given for cell twist distributions to 
which no immediate physical significance can he attached. These solutions 
are given in a brief manner to illustrate the mathematical methods that 
may he used to determine shear flows in terms of cell twistB . Distributions 
are chosep which appear to have the most promise of being useful in the 
development of a general theory of multicell sandwich plates - Since the 
twist function for a homogeneous, isotropic, elastic plate has a 
sinusoidal variation across the plate in certain natural buckle modes, 
a sine distribution of 0 may be of some usefulness in the analysis 
of multicell sandwich sections. 


Consider a sandwich section numbered as shown in figure 3( a ) • 
Assume the cell- twist distribution to be given by the follow in g formula: 


9 = a sin — . nx — r 

x 2 (el + l) 


(31) 


The general solution of the difference equation (equation (12)) becomes 
q x = B sinh 0x + C cosh 0x + b sin — - -y (32) 


where 


b = 


2AGa . 

1 - 2d cos — r ? r 

2(n + 1) 


] 


Za 


(33) 


and 0 is defined by equation ( 17 ) • 
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When the condition of antisymmetry and the "boundary condition are introduced, 
the solution "becomes 


•3x 



sin 


■six 

2(n + 1) 


slnh ftz 
sinh p(n + l) 


(34) 


The fomula has "been evaluated for a seven-cell member with square 
cells and equal wall thicknesses . Graphs of the cell-twist and shear-flow 
distributions are shown in figures 7(b) and 7(c)* A dashed line is drawn 
through the graphs at the centers of the cells to indicate the general 
distribution of the functions over the section. These continuous distri- 
butions are superimposed In figure 7(d) to indicate the difference in the 
distributions . The magnitudes of the cell twists and the shear flows are 
as follows: 


0 O 

0 1 

e 2 

0 3 

0 

0.383a 

0 .707a 

0.924a 


q 0 

q l 

*32 

q 3 

0 

0 .365b 

0.635b 

0.653b 


CHARACTERISTIC DISTRIBUTION OF 6 


When the static load distribution acting on a structure causes 
displacements such that the load and displacement are proportional at 
every point, the loads and displacements may be called characteristic 
loads and characteristic displacements • This designation was introduced 
in reference 4. The concept and designation may be applied to internal 
or external forces and displacements of a general nature. Arbitrary 
distributions of functions may oftentimes be conveniently esqpressed in 
terms of a series of characteristic distributions.. Hence, it becomes 
desirable to establish the characteristic-displacement functions that 
may be associated with a given type of structure. In the present case, 
characteristic distributions of 0 X and q x may be found by assuming 

0 x = ^ z (35) 

The difference equation (equation. ( 12) ) becomes 

■^Vl + (Sa ’ 2AGX)q x ' °2Vl = ° 


(36) 
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The solution of this homogeneous equation is given by 

q s = a sin 72; + b cos yx (37) 

The solution must be assumed in terms of circular, rather than hyperbolic 
functions in order to be able to satisfy the boundary conditions. Substi 
tution of this solution into equation (36) gives the following relation 
between 7 and X: 

X = fca - 2agCOB 7) 

= §5 (1 “ 2d 008 7) ( 38 ) 

Either term in the solution as given by equation ( 37 ) can be made 
to satisfy the boundary conditions by an appropriate choice of the coordi 
nate system. It will be assumed that the cells are numbered frcm 1 to n 
from left to right, as shown in figure 3 (c) • The first term can be made 
to satisfy the boundary conditions alone by a correct choice of 7. 

Hence, the coefficient of the second term may be set equal to zero to 
give 


b = 0 


and 


q x = a sin yx 


The boundary condition gives 

qn +1 = a Bln 7 (n + 1) = 0 


( 39 ) 

( 40 ) 


All functions connected with the problem are completely defined 
by n numbers. Hence the functions are said to exist in an n- point 
space . In such a space there Ib nothing of a useful nature to be gained 
by defining more than n distinct values of characteristic numbers 
such as 7 or X • There are n distinct values of 7 less than jt 
which will bring about compliance with the boundary conditions as 
stated by equation ( 40 ) . The distinct values, of 7 and X may be 
indicated as 7^ and X^. They are given by 


tot 

(k = 1, 2, 3, ... n) 

( 41 ) 

n + 1 

— (l 

- 2d cos — — — ) 

( 42 ) 

2 AG l 

n + 1) 
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and the resulting distribution 

(k = 1, 2, . . . n) (43) 

(k = 1, 2 ; . . • n) (44) 

* "k - - 

In equation (44) the ratio &kAk has taken the place of the coefficient 
in equation (39) • 


ARBITRARY DISTRIBUTION OF 0 


The n formulas represented by equation (43) or equation (44) define 
a closed orthogonal set of functions for an n-point space • Any arbitrary 
function may he represented by a finite Fourier series using these 
functions. If 0 X is an arbitrary distribution of the cell twists. It 
may be expressed as 


A characteristic distribution of 0 
of q are thus given by 


9 X = a^sin 
a k 


a = 




krtx 
n + 1 

krtx 


= 2^ a k 0 in (x = 1 , 2, . . . n) (45) 

k=l. 

The Fourier coefficients a^ remain to be determined. From equation (44) 
it is seen that the solution for shear flows is given by 

(x = 1, 2, . . . n) (46) 


This equation gives immediately the shear-flow distribution if the 
coefficients a;k for the cell twists are known. Equations (45) and (46) 
illustrate the fact that, in general, 6 and q are not proportional ■ 
functions . Although the individual terms of the two series are propor- 
tional, the sums of two such series are, in general, not proportional. 

The formula for the Fourier coefficients may be conveniently 
developed by using matrix algebra. Equation (45) represents a system 
of n equations defining a linear relationship between the 0^ values 

and the a ^ values . In matrix form this system is written as 

[0j = [S] [a] (47) 

In this equation [©] and [aj are column vectors (column matrices) 


<3x - 


i sin 


tatx 
n + 1 
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containing 0 X and a^. 'values . The square matrix [sj conta ins sine 
functions defined "by 

H * [°ik] 

^ m 

Consideration of the elements of [sj showB it to "be a symmetrical 
matrix having columns which are orthogonal vectors, each having a norm 

of Consequently, when H is multiplied hy itself, a scalar 

matrix is obtained. Thus, 

[ S ] [s] = [13 (49) 

where [I] is 'the identity matrix. The validity of equation (49) depends 
upon the orthogonality of the .columns of Js] . This orthogonality is 
expressed in scalar form hy the equation 


sin — ~ sin = S -f — 1 

4 n -+ 1 n + 1 2 


(i = J) 


) sin JS. sin -te- = 0 (ij*j) 

/ , n + 1 n + 1 
x=l 


This formula corresponds to the better known formula for the orthogonality 
of continuous functions . 


Multiplying through equation (47) by [s] gives 

[s] [ 9 ] - [s] [si [a] 


fi-PH 


M - H 0 


The matrix equation may be eapressed in scalar fozm to give 


„ , kjtx 

0-rSin — — - 
x n + 1 
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For any distribution of 0 X equation ( 53 ) gives the corresponding 
Fourier coefficients- These coefficients are then substituted into 
equation (46) to obtain the solution for shear flows- 


Since the shear flows are related linearly to the Fourier coefficients, 
which in turn are related linearly to the cell twists, there must "be a 
linear relation between the shear flows and cell twists . This relation 
may be derived and, thus, the necessity of computing Fourier coefficients 
may be avoided- The development of the relation may be most conveniently 
perforated with matrix algebra. 


Equations (45.) and (46) define two systems of equations between 
which the a^. values can be eliminated to obtain a single system 
relating 0 X and q x - Equation (45) bas been expressed in matrix form 
as equation (47) • In order to write equation (46) in matrix form, it is 
necessary to define the diagonal matrix [A] as 



X 1 

0 ... 

0 

[A] = 

0 

• • • 

0 


0 

0 ... 





m 


— 





m 


Equation (46) may be written in the following matrix form: 

[q] = TS] [A ]" 1 [a] 


(55) 


The reciprocal of a diagonal matrix is obtained by merely inverting the 
individual diagonal elements. 

Between equations (47) and ( 55 ) the column vector [a] may be 
easily eliminated - Solving equation (47) for [a] gives equation (52). 
Substituting into equation ( 55 ) gives , 

P -! ' 

[q] = [a] {s] [0] (56) 


For a given number of cells the matrix [S] is a standard universal 
matrix while the structural dimensions enter into the diagonal matrix [A]. 
The relative distribution of shear flows is dependent only on the relative 
aspect ratio of an internal wall segment as defined by the factor d- 
The solution given in equation (56) may be regarded as a general solution 
of the system of difference equations that govern the shear flows. 

Further discussion of a similar difference problem from this viewpoint 
has been given in reference 5* 
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As a numerical example , consider a section of four square cells with 
all vail thicknesses equal- The factor d equals O.25 and n equals 4- 
The matrices [S] and [a] becar.-o 


[S] = 


0.5878 

0.9511 

0 .9511 

0.5878 

0.9511 

0.5878 

-0.5878. 

-0.9511 

0.9511 

-0.5878 

-0.5878 

0.9511 

0.5878 

-0.9511 

0.9511 

-0.5878 


[A] = 


la 

2AG- 


0 .5955 
0 
0 
0 


0 

0.8455 

0 

0 


0 

0 

1.1545 

0 


0 

0 

0 

1.404-5 


Substitution of these matrices into equation (56) gives 

2.6796 O.7178 0.1914 0.0479' 
O.7178 2.8711 O.7657 0.1914 
0.1914 0.7657 2.8711 0.7178 
0-0479 . 0.1914 0.7178 2.6796 


V 


<*2 

q 3 

_ 4AG 
5Za 




e. 


e 


3 

0 4 


(57) 


(58) 


(59) 


For any given cell twists the shear flows may be immediately computed- 


It is possible to reduce the matrices to approximately half size in 
analyzing for arbitrary pell- twist distribution by taking advantage of 
the symmetry of the structure. The twist distribution may be divided 
into symmetrical and antisymmeirical parts. Separate matrix formulas 
such as equation (5 6), may be established for each part. The symmetrical 
and antisyrametrical parts of the shear- flow distribution would then be 
determined separately - 


T0RSI0N-C0HSTAHT F0KM0IA 


By the use of the formulas that have been previously developed for 
shear flows in pure torsion, it is possible to derive a formula for the 
torsion constant for a multicell sandwich section. . (The author is indebted 
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to M. A. Biot for assistance in the development of this formula.) The 
torsion constant may he introduced hy expressing the relationship "between 
section torque of a beam and rate of twist at that section as 

T = GJS (60) 

where G is the shear modulus of elasticity and J is the torsion 
constant. Consider a section with an odd number of cells numbered 
symmetrically as shown in figure 3(a)* The total torque acting on a 
section is eq ual to twice the volume beneath the stress surface and hence 
is given by 


T = 2A ^ g x 


Ls 

-n 


(61) 


The formula for shear flow, as given by equation (21 ) , may be substi- 
tuted into equation (6l)j thus. 


T = 2AD 


\ i-, cosh gx 1 

_Z_. L cosh 0 (n + 1 ) J 


= 2(2n + l)AD - 


2AD V 

cosh p(n + l) Is 
-n 


cosh 0x 


(62) 


The summation contained in the second term. 'must be evaluated. For this 
purpose it may be separated as follows: 


-1 


n 


y cosh 0x = y cosh 0x + 1 + y cosh 0x 
Introduc in g exponential functions gives 1 


(63) 


cosh 


V (e Px + e^ 1 ) 

n 


(64a) 

(64b) 


- n 

^ cosh 0x 
1 1 

Formulas (64a) and (64b) may be regarded as finite geometric power 
series in e^ or e - ^. The following known formula may be used: 


a 2 = 


a(l - a 11 ) 


1 


1 - a 


(65) 
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From this formula it is seen that 

n. -1 

^ cosh Px = ^Tcosh px 
1 -n 


e^(l - e n P) + e~Kl ~ e~ n ft) 
1 - e P 1 - e"P 


Substitution in eguation (63) gives 


n 

I 


cosh Px 


= 1 + eP(.l - e 11 ^) + e~P(,l - e~ D ft) 


1 - eP 


1 - e“P 


This formula may be condensed as follows: 
n 


V cosh PX = 1 + - ^5+ e^ / 2 ^ 1 - e-n pj 

L e -P/2 _ e P/2 e P72 _ 0 -P/2 


,H)P ■ e~(4> 

e P/2 : e'-p'/s 


• sinh (2n + 1)| 

sinh ■& 

2 

This formula may now be substituted into equation (62) to give 

2AD sinh ( 2n + l)£ 

T = 2(2n + l)AD 

cosh (n + l)p sinh q 

= 2(2n + l)AX^(p,n) 

where 

sinh (2n + l)£ 

<f = 1 5 

(2n + l) cosh (n + l)p sinh £ 

For calculation purposes it is somewhat more convenient to express 
the follow i ng equivalent form: 


(66) 

(67) 

( 68 ) 

( 69 ) 

(70) 
$ in 
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2 sinh | 

r? 

♦*> 


(2n + l) 

■a 

•ri 

1 m 1 

(f + U 

Ip - sinh | 

( n + i>] 


0 = 1 " 


Substitution of the formula for D from, equation (l6) into 
equation (69) gives 


(71) 


m _ 4(2n + l)A 2 G60(P.n) 

(l - 2d)Za 

This formula may be substituted into equation (60) . Dividing out the 
factors G0 gives the following formula for the torBion constant: 


T 4(2n + l)A g 0(P,n) 
J " (1 - 2d)Za 


(73) 


It is convenient to introduce the tors ion- constant formula for a single 
cell as given by 


J 1 - ^ <74-) 

If K is the number of cells (E = 2n + l) , the torsion constant for a 
section having E cells is given by 


E 



0(P,n) " 
1 - 2d 


(75) 


Since P is determined when d is known, as shown by equation (17), 
0 may be regarded as a function of d and E rather than p and n. 
Thus, it becomes convenient to introduce 7|r(d,E) as 


ty(d E) = -^(P-i. n ) . '0(l»H) 
’ 1 - 2d 1 - 2d 


(76) 


The torsion constant becomes 

J E = EJ 1 ^(d,E) 


The function ^ is shown graphically in figure 8 for a section 
having d = 0 .25 . 


(77) 
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CONCLUSIONS 


Fonmulas have 1)0011 developed in explicit form for the ahear flows 
in multicell sandwich "beams in "bending and torsion- Formulas have also 
"been developed for the shear flows corresponding to various functional 
chordwise distributions of the cell twists and, in scalar or matrix 
form, for the shear flows corresponding to an arbitrary chordwise twist 
distribution. A formula for the torsion constant of a multicell sand- 
wich beam has also been developed. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., September 9, 1948 
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(a) Odd number of cells with 
Symmetrical numbering. 



(b) Even number of cells with 
' symmetrical numbering. 



(c) Arbitrary number of cells with 
unsymmetrical numbering. 


Figure 3.“ Cell numbering. 
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/ 


Figure 5.- Imaginary cell used 


in statement of 
boundary oondi+ion . 



(b) Cell numbering, 



Cc) Shear -flows. 


Figure 6.- Shear -flows due +o bending 


















